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Abstract 

The fluid models mentioned in the title are classifled. All char- 
acteristics of the fluid are expressed through a master potential, sat- 
isfying an ordinary second order differential equation. Different con- 
straints are imposed on this core of relations, flnding new solutions and 
deriving the classical results for perfect fluids and dust as particular 
cases. Many uncharged and charged anisotropic solutions, all confor- 
mally flat and some uniform density solutions are found. A number 
of solutions with linear equation among the two pressures are derived, 
including the case of vanishing tangential pressure. 

1 Introduction 



The description of gravitational collapse and evolution of compact objects 
under various conditions remain among the important problems of general 
relativity. They are described by spherically symmetric relativistic fluid mod- 
els where the metric depends on time and radius. In general, they possess 
shear, expansion and acceleration, which makes them rather hard to study. 
In some physical situations, however, perfect fluid (PF) models without accel- 
eration or dust are a good approximation. The general dust solution, known 
also as the LTB solution [Ij , ^ , [3] , |1] is one of the most exploited examples. 
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It has three branches, which are characteristic for all other solutions. Much 
later the PF solution for the parabolic branch was found [5], [B]. It was shown 
[7] that the other two branches lead to the Emden-Fowler equation [8] and 
some exact solutions were proposed [9] . Solutions which obey an equation of 
state have constant pressure and density, which is equivalent to the presence 
of cosmological constant [TU]. The case of charged dust has also attracted 
attention jJT], [12], [IS]- Such solutions cannot be shear-free [T^. Finally, 
the general solution was found in mass-area coordinates [T^. Recently, a 
charged geodesic PF solution with constant pressure was studied [15|. The 
LTB model has been generalized to include dissipative fluxes [16] . 

On the other hand, different mechanisms have been identified through the 
years which create pressure anisotropy in stellar models and make the fluid 
imperfect [17]. Such are the exotic phase transitions during gravitational 
collapse, the existence of a solid core or the presence of a type P superfiuid. 
Viscosity may also be a source of local anisotropy as well as the slow rotation 
of a fluid. It has been shown that the sum of two PFs, two null fluids or 
a perfect and a null fluid may be represented by effective anisotropic fluid 
models [IS]. Recently it was pointed out that the same is true for PF with 
charge, bulk and shear viscosity [19\. These arguments have stimulated the 
study of anisotropic fluids. Geodesic anisotropic spheres were discussed in 
non-comoving coordinates and a particular solution was given [20]. Lately 
the notion of Euclidean stars was introduced and their properties were inves- 
tigated with or without radiation [H], [22] • In the first case the fluid motion 
is geodesic. The evolution of the shear in geodesic fluids was studied too 
|23j . Unfortunately, the extensive classification of Einstein solutions [4] in- 
cludes only uncharged PF. Static charged PF solutions have been described 
elsewhere [21] 

The aim of the present paper is to generalize the above results and give 
a classification of the geodesic anisotropic fluid solutions, similar to that 
of shear- free anisotropic fluid solutions [25]. For this purpose we use the 
reformulation of the Einstein equations in terms of the mass function [I], 
[2S] • In this approach, besides the new solutions, the classical results reviewed 
above follow as particular cases. 

In Sec. 2 the field equations based on the mass function are given for 
non-radiating spheres in geodesic flow. The fundamental potential and the 
main equation that it satisfies are introduced. Expressions for the other 
characteristics of the metric are given. Sec. 3 discusses the general and some 
concrete solutions of the geodesic anisotropic model. In Sec. 4 the particular 
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cases of PF and dust are investigated, making contact with the existing 
hterature. In Sec. 5 the general charged solution and its subcases of charged 
PF and charged dust are discussed. Sec. 6 studies the case when the radial and 
the tangential pressures are proportional. In Sec. 7 the general conformally 
fiat solution for geodesic anisotropic fluids is found. The uniform density 
solution is studied in Sec. 8. In Sec. 9 a comparison is made between shear- 
free and geodesic fluids. The final section contains some discussion. 



2 Field equations for geodesic fluids 

Spherically symmetric relativistic fluid models are described by the metric 

ds^ = -e^-de + e^Hr^ + {dd^ + sin^ dd^^) , (1) 

where u, A and R are functions of t and r only. The spherical coordinates are 
numbered as x° = t,x^ = r, = 6,x^ = (p. The Einstein equations read 

- \9^JR^ = SvrT,,. (2) 

Here Rij is the Ricci tensor, Tij is the energy-momentum tensor and we have 
set G = c = 1. For an anisotropic fluid model without radiation one has [27] 

Tik = (p + Pt) UiUk - Pt9ik + (Pr - Pt) XiXk- (3) 

Here p is the energy density, Pr is the radial pressure, pt is the tangential 
pressure, is the four-velocity of the fluid, is a unit spacelike vector 
along the radial direction. These vectors satisfy the relations 

u'ui = -1, x'Xi = 1, u'x^ = 0. (4) 

The coordinates are assumed to be comoving, hence 

u' = e-'^Sl x' = e-%. (5) 

The usual field equations (2) are rather cumbersome except for the (01) 
component 

R' - Ru' - R'X = 0. (6) 

The dot above means time derivative, while the prime denotes a radial one. 
It is easier to work in the mass function formalism [1], [26], namely 

m' = AnpR^R', (7) 
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rh = -A-aprR R, (8) 

Hp + Pr) = -P-^{P + Pt), (9) 
^' {P + Pr) = -P'r + ^ {Pt - Pr) . (10) 



The mass function is given by 

1 



m = -R i^e-'^R' + kj , (11) 

k = l- e-^^R'\ (12) 

where k is the binding energy [28j. Bound configurations have < < 1, 

unbound ones have k < and the marginally bound case is k = 0. These 
three cases are called also eUiptic, hyperbolic and parabolic respectively. The 
second Weyl invariant reads [25j 

/ m \ ' R 4:71 ^ , , . 
The expansion of the fiuid is 





e = e-^^AH--J, (14) 

while the components of the shear tensor are proportional to 

(15) 

The four acceleration is iii = (0, z/', 0, 0). 

Geodesic fiuids have vanishing acceleration, hence z/' = and i/ is a 
function of time only. It can be set to zero by a coordinate transformation. 
If = Eqs (6) and (12) lead to geodesic fiow. The contrary is not true. 
When u = and R' the same equations yield k = k{r) and 

E)/2 
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It is not possible to make k vanish by a transformation of r. The so called 
Euclidean star models pi] have k = 0. Eq (8) becomes with the help of Eq 

(11) 

(^RR^^ ' = -MprR^R - kR. (17) 

This is the main equation for geodesic fluids. It is second order and contains 
only time derivatives. The dependence on r is viewed upon as parametric. 
One gets rid of the first derivatives by the substitution 

R = Z^/\ (18) 

which results in 

Z = -6nprZ - ^kZ-^/\ (19) 

We can find Z when pr is given or vice versa. The second case is much simpler 
but it is not guaranteed that the resulting p^ will be physically realistic. Thus 
the general solution depends on one arbitrary function of t and r. 
Eq (10) becomes 

p'r = ^iPt-Pr) (20) 

and gives an expression for the tangential pressure 

P. = (21) 

When Pr = 0, pt vanishes too and we get a dust solution, but the opposite is 
not true. The mass is found from Eq (11), which becomes 

m = ^R + A;) = ^Z^ + ^Z^/^. (22) 

Then the energy density follows from Eq (7) 

4(zA' + 9{kZ^/'y 
24vrp = ^^-J^' • (23) 

The expansion and the shear scalar are given by 

e^ff. (24) 
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3(lnZ)' 

A geodesic solution is also shear-free when Z and consequently R is separable. 
It has no expansion when R is of the form 

i^(^,r) = [i^l(^) + ^2(r)]'/^ (26) 



3 General and particular solutions 

The main equation (19) may be simplified further by the transformation 

Z = iiW, ^ = /i^ (27) 

where /z (t, r) determines the new time variable r. It becomes 

W^r = - (/i + 67rp,/i) /iW - ^A://3yf/-i/3^ ^28) 

We choose /i so that the bracket vanishes and the above equation breaks into 
two 

6vrp. = = ^/iz^, (29) 
/i jjp 

W,, = -^A;/i«/3iy-V3. (30) 

The general solution follows when one chooses arbitrary k{r), n{t,r) 
or pr{t,r), and finds W{T,r) from the above equation. This determines 
consecutively Z {t,r), R(t,r) and all other characteristics of the fluid. Eq 
(30) for k ^ resembles the Emden-Fowler (EF) equation [8] 

Wrr = At^'W^ (31) 

with A = —jk and / = —1/3. In order to obtain analytical solutions we set 
/i = T^"'/^. Eq (27) yields for n ^ —4/3 and n = —4/3 respectively 

3n+4 3n + 4 , , , „ 

r— = ^^(t-a), r = e*-", (32) 

where a (r) is an arbitrary function of integration. Eq (29) gives for n ^ —4/3 
and n = —4/3 respectively 



There are general solutions of Eq (30) when n = (/i = 1) or n = —8/3 
(/i = 1/r) [H], but according to Eq (29) they are particular dust solutions, 
which are discussed in the next section. The case n = —4/3 (/i = r~^/^) 
is also soluble, but leads to constant pressures, representing effectively a 
cosmological constant. 

There is, however, a particular two-parameter solution of Eq (31) for any 

(n + 2)(n + / + l)" 



n + 2 

W = ar -^-i a 



A{l-lf 

In our case / is fixed, while n may be a function of r. We obtain 



(34) 



/ N 9 (3n + 4)^/c 

R = ai(t-a), a? = H rr^ r- (35) 

^ ^ 4(n + 2)(3n + 2) ^ ^ 

The solution depends on three arbitrary functions of the radial coordinate 
k, a, n. When a = 0, i? is separable, hence, the solution is shear-free in 
addition to being geodesic. 

No other explicit solutions exist when is a power of r. However, there 
is one with /i = (ar^ + 6r + c) where a, 6, c are functions of r . It leads 
to time- independent radial pressure pr {f)- The pressure profile is given by 
the initial conditions and may be arbitrary. This case is better approached 
from Eq (19), which is integrated simply by multiplication with Z. We get 
the integral 

dZ , , 

^^^=^^^= = t + s, (36) 

-QnprZ^ - \kZ^I^ + h 

where h{r) ,s (r) are arbitrary integration functions. The first is related to 
the mass functions through Eq (22) 

4 2 

m = --nprR^ + -h. (37) 

The integral in Eq (36) has an analytical expression when any of the functions 
Pri k, h vanishes. Thus, when h = we resort back to R to find 

R={- sin-v/6^(t + s), (38) 

ySirprJ 3 
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provided Pr > and A; < 0. When A; = we obtain in a similar way 



^ h^^^^y^^ (39) 

This is a concrete example of an Euclidean star. Finally, when = a dust 
solution follows, to be discussed in the next section. 

Up to now we have assumed that A; 7^ in Eq (30). The case k = gives 
immediately 

W = aT + b, T= fi^'^dt, (40) 
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where a{r) ,b (r) and fi {t, r) are arbitrary. The radial pressure is found from 
Eq (29). When b = 0, R becomes separable and the solution is shear-free 
too. 



4 Perfect fluid and dust 

Perfect fluid may be considered as anisotropic fluid with the equation of state 
(isotropy condition) = = P- In this case Eq (20) shows that p = pit). 
It is enough to take a /i (t) in Eq (29). In fact, the transformation given by 
Eq (27) was applied first to PF [7], [9]. The solution (36) was also proposed 
[9], without noticing that it holds for p (r) and contradicts the PF condition. 
The concrete solution for A; 7^ described by Eqs (33, 35) holds for PF too, 
provided a,n are constant. The general solution in the case A; = is given 
again by Eq (41), but with /i uniform in space [5], [6]. Szafron [29] has found 
the particular solution /i = with q any real number. 

When the pressures vanish, one has collapsing dust. Eq (10) shows that 
its motion is geodesic (without acceleration). According to Eq (8) the mass 
function is time- independent. Its profile characterizes the solution. The 
main equation (19) may be integrated once and becomes Eq (22). A second 
integration leads to Eq (36) with = and h = 9m/2 , in accord with Eq 
(37). This is the well known LTB solution [1], [2], [3]. Eq (22) coincides with 
Eq (73) from Ref. [16], where one can find the analytic expressions for R 
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in the parabolic, hyperbolic and elliptic case. For example, when k = one 
easily finds 

+ ' ^^^^ 
m{r) ,h (r) being arbitrary functions. The solution is also shear-free when 
6 = 0. 



5 Charged fluid 



Spherical symmetry allows the appearance of only an electric field E in the 
radial direction. The energy-momentum tensor of this field can be described 
[in] as addition of the following effective energy density and pressures to the 
fluid 

P^=pf = -Pf = f. (43) 
The Maxwell equations give 

- = 47rC7 = Ee-\ (44) 

where C is the charge function of the fluid and q (r) is an arbitrary function. 
The main equation in its two forms (19) and (30) becomes 

Z = -GTTprZ - ^kZ~^/^ + 6nqZ~^/^, (45) 

Wrr = -^kfi'/'W-'/^ + 67rqfi^/^W-'/\ (46) 
while Eqs (27) and (29) do not change. Eq (21) transforms into 

.=i^-^ U7) 
^* (i?2)' 2R^R' ^ ' 

Now when p^. vanishes pt is not obliged to vanish, giving a class of non-dust 
solutions. 

For 7^ Eq (46) becomes the modifled EF equation 

Wrr = Ait'^'W'^' + A2r"W^2 (48) 
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when /i = r^"^/^. Then 

15 3 

mi = --, ma = --, Ui = 2n2, Ai = --k, A2 = Qirq. (49) 

There is a table of 108 exphcit solutions of Eq (48) [8]. Those which 
satisfy the conditions above have rii = 0; — |; — |. The first two lead to 
vanishing Pr, the third has it constant. They are non-trivial and non-dust. 
They can be approached also as particular cases of the class of solutions with 
Pr (r). Like for uncharged fluid one integrates directly Eq (45) to find 

3 f RdR , . , ^. 

^/ I ^=t + s{r). (50) 

^ -J ■ l-Q-nprR"^ - \kR^ + hR- lOvrg 



Eq (50) coincides with Eq (36) if g = 0. When h or pr vanishes, the integral 
is expressed in terms of elementary functions in five different ways depending 
on the relations between k, q and Pr or h. 

For k = 0, Ai = and the first term in Eqs (46, 48) disappears. We get 
once more the usual EF equation hke Eqs (30, 31) in the uncharged case, but 
the power of W is different. This change brings a host of explicit solutions 
with n2 = 0; — |; — |; — y ; — |; — |; — |; 1; 2. The first two have vanishing 
radial pressure. The particular solution (34) also holds provided / = — 5/3 . 

Charged PF has pr = Pt = P- Eq (47) becomes 

This is another relation between p and R in addition to Eq (45). A possible 
solution is p (t) (like in the uncharged case) and q = const. Then n = nit) 
and a, n should be constant in the above solutions. The isotropic pressure is 
given by Eq (33). 

When p' 7^ one should replace R from Eq (51) into Eq (46) having in 
mind that 

W = li-^R'/^ (52) 

Even in the simplest case /c = this leads to contradiction, because the result 
is n = —4/3 and hence, p = const. 

Charged dust may be considered as a subclass of the PF solutions with p = 
0. The general solution is known, though in another coordinate system 
However, Eq (10) does not lead now to u' = because of the electromagnetic 
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additions to the energy density and the pressures given by Eq (43). We have 
to impose this condition in order to make the fluid flow geodesic and the 
solutions found here form a subclass of the general Ori solution. When p = 
Eq (51) yields q (r) = qo = const. It is not necessary to introduce W and /x. 
Eq (45) becomes 

Z = -^kZ-^^'' + 67rgo^"'/'. (53) 

This is the modified EF equation (48, 49) with ni = n2 = and A2 = Gnqo- 
Its solution is given by example 2.4.2.56 from [8] and has three branches. 
When k <0 

t = Cie"^ + Cae-"^ + C3A, R = u (Cie"^ - Csg-"^) + C3, (54) 



where A is a parameter, u = y and the integration functions Ci (r) satisfy 
the relation 

3 (AiC| + A2) + I6AIC1C2 = 0. (55) 

When k > 

t = Ci sin + C2 coswA + C3A, i? = w (Ci coswA — C2 sincuA) +C3, (56) 

3 + A2) + AA\ {Cl + Cl) = 0. (57) 

Finally, when /c = Eq (53) may be integrated by multiplying with Z. 
The result is 

^'^^ lt + C2, (58) 



Y 
where 

Y = {CiR-l^Tiq^f^ . (59) 
Eq (58) becomes a cubic equation for Y 

Y^ + bAT^q^Y = Cl [t + '^-C^ . (60) 

It may be solved for Y and R. The solutions are not separable and therefore 
are not shear- free in accord with Ref. Illl. 
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6 Solutions with pt = ^Pr 



This equation of state among the pressures includes the cases of perfect fluid 
and Pt = 0. Eq (21), which in principle determines the tangential pressure, 
now gives an expression for pr in terms of Z 

Pr = h (t) zl^^-'\ (61) 

where h is an arbitrary positive function. The main equation (19) becomes 

Z = -h{r)Z-'/^ -QTTh{t)Z^. (62) 

This represents the modified EF equation (48) with 

47 - 1 3 
nil = -1/3, m2 = — - — , rii = 0, Ai = --/c, A2 = -Gn, (63) 

provided h = t"-^. The following explicit solutions exist: 

When 7^ we may take ^2 = 0, m2 = —5/3 {h = ho,^ = —1); then Eq 
(62) coincides in form with Eq (53) and the solution (54-57) holds. Other 
possibilities with the same m2 are n2 = 1 and n2 = 2 which lead to examples 
2.4.2.58 and 2.4.2.54 respectively [8]. The case pt = (7 = 0) reduces Eq 
(62) to 

Z = -(^k + 67ihj Z'^/\ (64) 

It is easily integrable for constant h = ho when it coincides in form with the 
LTB solution. 

When k = Eq (62) becomes the usual EF equation with mi = ^-^^ 
provided h = t"^ One can choose from the 28 analytic solutions (5 one- 
parameter families and 23 isolated points) |S] and the particular solution 
(34). 

Eqs (21, 61) show that when 7 < 1 (including the case of vanishing 
tangential pressure) the radial pressure becomes infinite at i? = 0, adding 
another singularity to the curvature singularity of the LTB solution. This is 
true provided there is no bounce of the collapsing fluid and the point R = 
is covered by it after some time. A similar situation is described in Ref. [28] . 
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7 Conformally flat solutions 

These solutions have ^'2 = 0. Combining Eqs (13, 20) yields for geodesic 
fluids 

(13+2-^31^ (65) 

The vanishing of the second Weyl invariant requires 

TTl 

— + 27rpr^f{t), (66) 

where / is an arbitrary positive function since m and Pr are positive for 
realistic star models. Eq (8) transforms into 

^ ^ X^f-2npr. (67) 



2f + X R' 
Inserting Eq (22) into Eq (66) results in an expression for p^ 

2vrp. = /-^^. (68) 
We combine this equation with Eq (66) to find 

2RR -if^ -^fR^ + k, (69) 

which is analogous to Eq (17). The first derivative is eliminated by the 
substitution R = Y'^ 

Y = -fV + (70) 

This is the main equation in the conformally flat case, replacing Eq (19). It 
has no reference to the pressure, although / enters it in a similar way. Once 
again we use the transformation Y = /iW which splits Eq (70) into two parts 

/i + /// = 0, Wrr^^W-^ (71) 

where fi = fi (t). Once again we obtain a EF equation. This time it is easily 
integrated by multiplying with Wr- We get 

R = bf,' + V^J fi-'dt + c. (72) 
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Here b{r) ,c (r) are arbitrary integration functions and /i^ is an arbitrary 
positive function of time. This formula is valid for any k. In the marginal 
case / = (which requires some radial tension) it simplifies to 



(73) 



The general conformally flat geodesic solution is also shear-free provided 

c = = A;. 

Let us find now the conformally flat PF solution. Eq (20) shows that 
p = pit). Eq (68) becomes 



X{t) 



Eq (22) then gives for the mass function 

m^^X (t) R^ 



(74) 



(75) 



and Eq (7) shows that the energy density is uniform, p = p{t). Putting the 
general solution (72) into Eq (74) yields 



1 + 2pfiV 

1t/2 I k 



X, (t) . 



(76) 



The l.h.s. should depend only on time, which is ensured by setting to con- 
stants c = Cq and k/b"^ = Ci. R becomes 



R^bjj? 



H dt-\- Cq] + Ci 



(77) 



The solution is shear-free because R is separable. 

The conformally fiat dust solution is a subcase of the PF solution with 
vanishing pressure. Eqs (67, 68) become 



3/ 



R 

R' 



They lead to 



R'^ + k^ 2fR^. 
R^g{r)f{t)-"\ 



(78) 
(79) 
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/ = 3^/^2/1/3-^, (80) 

where g (r) is arbitrary, but positive. Obviously k/g"^ = C2 = const. The 
equation is integrated by the substitution y = f^/^ and the final result is 

/ o \ 2/3 

R = [^) 9 ir) (t + cf' , (81) 

where C3 is some constant. This formula resembles the k = case of the LTB 
solution given by Eq (42). Like the PF solution (77), the dust solution is 
shear- free. 



8 Uniform density 



The energy density in this case is a function of time, p (t). Eq (7) gives upon 
integration 

3 ^ ' 



(82) 



where an integration function has been set to zero |30j. Plugging this formula 
in Eq (22) yields 

Stt R' + k 
With the help of Eq (8) it is transformed into 



3k = R^ 



iirp — 



(84) 



S{p + Prf. 

If we put Eq (82) into Eq (8) an expression for the radial pressure is obtained 



R Air 

- '"^^ = ^ + Y^- 



(85) 



Let us discuss first the case k = 0. Eq (84) defines Pr in terms of the 
arbitrary p (t) 

Pr = - P- (86) 



V24^ 
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It shows that Pr — Pr (t) and the fluid is perfect. Eq (83) deflnes R 

R^a{r) exp j ^rpdt^ (87) 

with integration function a (r). The solution is shear- free. 

In the dust subcase Eq (86) becomes an equation for p. Its solution is 

^ ~ 67r(t + to)" ^^^^ 
where Iq is an integration constant. Then R simplifles 

i? = a (r) (t + Uf'^ . (89) 

It coincides in form with the conformally flat dust solution given by Eq (81). 
However the function g (r) there depends on k and it holds for any k. 
Let us consider now the case k 0. We set 

R {t, r) = ^/\k{r■) \P {t, r) . (90) 

Eq (83) becomes 

This equation is not integrable in general and separable solutions are not 
possible. For PF, however, Eqs (84, 90) show that P = P (t) and we may 
take it as an arbitrary function which determines p and R. The pressure is 
found from Eq (85). In the dust subcase m = m{r) and Eq (82) becomes 

The l.h.s. depends on t, while the r.h.s. depends on r, hence, both of them 
are constant, 3co/87r. Then Eq (91) reads 

= (93) 

and may be integrated, giving a rather complicated inexplicit expression for 
P 



'P (co - + 2 1^ 



P - ecQ 

It is not possible to recover the A; case by setting £ = 



e{t + to) . (94) 
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9 Comparison with shear-free anisotropic fluid 



The structure of the present classification of geodesic anisotropic fluid spheres 
is analogous to that of shear-free anisotropic fluid spheres [25], denoted fur- 
ther as I. In both cases the (01) component of the Einstein equations may 
be integrated. As a result all metric components depend on R plus arbitrary 
functions of one variable like G {t) in I and k (r) here. The latter function 
divides any solution into two or three branches, namely k = and k 
{k > and < in some cases). 

In both cases the fundamental ingredient is a second order differential 
equation for R in only one of the variables t or r. It is based on the expression 
for the mass function, which luckily contains derivatives of only one kind. For 
shear-free fluids this is Eq (I 34) with only first order derivatives in r. It leads 
to the main Eq (I 36), which is of second order. For geodesic fluids the mass 
function is given by Eq (22), containing only a time derivative. It leads to 
the main Eq (19), which is of second order in time derivatives. All other 
characteristics of the fluid are expressed through the solutions of the main 
equation 

Shear, expansion and acceleration are usually expressed through the met- 
ric. An essential role in the formalism plays a second formula for the vanishing 
characteristic, which holds for any fluid sphere. It elevates its dependence 
from metric components to sources (pressures and energy density). In I it is 
given by Eq (I 29) for the shear, whose derivation required the mass decom- 
position formula (I 15). Here this is Eq (10) for the four-acceleration, which 
enters the usual set of Einstein equations, based on the mass function. It is 
not necessary to use the mass decomposition formula except for the definition 
of ^2 given by Eq (13). 

The main equation contains no first derivatives and just one non-derivative 
term, when one passes from R to L = r / R in I or to Z = R^^"^ here. For 
geodesic fluids a second step is necessary, namely the transformation in Eq 
(27). The general solution of the main solution is based on one arbitrary 
function. In I this may be R, L, Z, Ap or \l/2- Here these are R, Z or pr- 
When R is chosen, solutions are found quite easily, however, they may be 
unphysical and the passage to PF is rather involved. Therefore we take Z or 
Ap as arbitrary functions in I or pr here and solve the main equation. Ana- 
lytic solutions of it are found by reduction to an EF equation. This happens 
when an arbitrary function is taken as a power of in I or r in the present 
paper. 
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The passage to PF is controlled by Ap in I (it enters directly Eq (I 36)) 
and by p[. = here, which makes the radial pressure dependent on time 
only. For geodesic fluids the additional case of dust exists. It is treated as a 
subcase of PF, when the pressure vanishes. 

Charged fluid has additions to the pressures and density which bring a 
second term to the main equation. It becomes the modified EF equation and 
possesses a number of one-parameter or isolated solutions. 

There is one easily integrable charged case when some arbitrary functions 
are set to constants for shear-free fluids or pr = Pr (r) for geodesic fluids. The 
corresponding integral (I 52) leads to the Weierstrass elliptic function in I and 
to Eq (50) here, which has analytic expression in some cases. In I an equation 
of state for PF emerges from this case. Here the necessary condition for 
this is constant pressure, which represents a cosmological constant in accord 
with Ref.[TO] and is rather trivial. The case of charged dust is completely 
integrable. 

The general conformally flat solution is found in both papers and is given 
by Eq (I 64) and Eq (72) respectively. The solution for PF follows when an 
arbitrary function depending on r is put to zero in I, or two such functions 
are made constant for geodesic fluids. Conformally flat dust is obtained by 
further simpliflcation. Conformally flat PF solutions intertwine with uniform 
density solutions. 

The general uniform density solution is found in I, while here this is done 
in the k = case for anisotropic fluid and in all cases for PF and dust. 

An important case is to impose a linear equation of state pt = jpr between 
the two pressures, which is a generalization of the isotropy condition 7 = 1, 
leading to PF. In both cases there are analytic solutions for R at some discrete 
values of 7, including the case of vanishing tangential pressure 7 = 0. 

A geodesic solution may be shear-free too. This happens when R is sep- 
arable. Many of the above solutions become separable when the integration 
functions are chosen appropriately. However, we were unable to flnd solu- 
tions that are both geodesic and expansion-free, i.e. satisfy Eq (26), except 
for the trivial case R{r). 

As a flnal remark, the EF equation or its modifled version appear fre- 
quently in both the shear-free and the acceleration-free case. Integrals which 
lead to the Weierstrass function are as frequent in I but here they do not 
appear at all. 



18 



10 Discussion 



Recently we have stressed the importance of anisotropic fluid models for 
astrophysics [19]. The fluid flow has three important characteristics - shear, 
acceleration and expansion. It is hard to obtain general solutions when all 
three are non-trivial. We have classified in another paper the shear-free 
anisotropic spheres [25]. In the present paper the same is done for geodesic 
flows, where the acceleration vanishes. 

Taking the fluid anisotropic means accepting the maximum freedom al- 
lowed by spherical symmetry. The general solution depends on an arbitrary 
function and a core of simple relations is obtained. Only after that we start 
to impose different constraints aiming the system at particular cases. 

Thus, PF is this core with imposed isotropy condition pr = pt, charged 
fluid is a neutral fluid with a special kind of anisotropy, conformal flatness 
means the constraint \l/2 = 0, uniform density constrains the energy density 
to a function of time, anisotropy may be studied in more detail by a linear 
equation among the pressures, which includes the cases where one of them 
vanishes. In this way the relations among the particular cases become much 
clearer due to the relations of all of them to the core. Time and again 
the second order Emden-Fowler equations appear, which sometimes may be 
integrated to first order and even to implicit integral formulas for the solution, 
like Eqs (36, 50, 58). In many cases they are expressed through elementary 
functions. 

Imposing one constraint fixes the arbitrary potential of the general model 
but up to functions of time or radius. Imposing a second constraint makes 
the system overdetermined, yet special solutions still exist when the arbitrary 
functions become constants. In total, we find many uncharged and charged 
anisotropic solutions, all conformally flat solutions, a large class of solutions 
with proportional or vanishing pressures and some uniform density solutions. 
The general geodesic dust solutions are found explicitly in all cases, building 
upon the LTB solution. 

The present classification is based on the mass function formalism, which 
is applicable only to metrics varying with time. Therefore the static case 
should be studied on the base of the usual Einstein equations (2) which 
simplify a lot in this case. 

Much work remains to be done. One should classify as a next step the 
expansion-free anisotropic spheres. Physically realistic solutions for collaps- 
ing star models should be distinguished from the above classes, cases and 
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subcases with numerous arbitrary functions and constants of integration. 
Junction conditions give additional constraints. It seems possible to general- 
ize the above constructions also to radiating anisotropic spheres, which have 
attracted a lot of interest recently. 
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